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Abstract 

In this paper, we introduce a new bivariate distribution we called it bivariate expo¬ 
nentiated modified Weibull extension distribution (BEMWE). The model introduced 
here is of Marshall-Olkin type. The marginals of the new bivariate distribution have 
exponentiated modified Weibull extension distribution which proposed by Sarhan et 
al.(2013). The joint probability density function and the joint cumulative distribu¬ 
tion function are in closed forms. Several properties of this distribution have been 
discussed.The maximum likelihood estimators of the parameters are derived. One real 
data set are analyzed using the new bivariate distribution, which show that the new 
bivariate distribution can be used quite effectively in fitting and analyzing real lifetime 
data. 

Key words: Joint probability density function, Conditional probability density func¬ 
tion, Maximum likelihood estimators, Fisher information matrix. 


1 Introduction 

Recently, Sarhan et al. (2013) lias defined a new four-parameter distribution referred to 
as exponentiated modified Weibull extension (EMWE) distribution. Sarhan et al. (2013) 
defined the (EMWE) distribution by exponentiating the new modified Weibull extension 
(MWE) distribution which discussed by Xie et al. (2002) as wos done for the exponentiated 
weibull (EW) distribution by Mudholkar et al. (1995). They observed that exponential dis¬ 
tribution, generalized exponential distribution (1999), Gompertz distribution (1824), gener¬ 
alized Gompertz (GG) distribution (2013), exponentiated Weibull (EW) distribution (1995), 
Weibull extension model of Chen (2000), modified Weibull extension (MWE) distribution 
(2002) and etc distribution can be obtained as special cases of the (EMWE) distribution. 

The objective of this paper is to provide a new bivariate distribution, whose marginals 
are (EMWE) distributions which referred to as bivariate exponentiated modified Weibull 
extension (BEMWE) distribution. It is obtained using a method similar to that used to 
obtain Marshall-Olkin bivariate exponential model Marshall and Olkin (1967). 
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University, Mansoura 35516, Egypt. E-mail address: elgoharyO@yahoo.com (A. El-Gohary) 
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The paper is organized as follows. Section 2 presents the shock model yielding the (BE- 
MWE) distribution. Also, the joint cumulative distribution function, the joint probability 
density function, the marginal probability density functions and the conditional probabil¬ 
ity density functions of (BEMWE) distribution is derived in Section 2. In Section 3 sum 
reliability studies are obtained. Section 4 presents the the marginal expectation of the (BE¬ 
MWE) distribution. Section 5 obtains the parameter estimation using MLE. In section 6 a 
numerical result are obtained using real data. Finally, a conclusion for the results is given 
in Section 7. 


2 Bivariate exponentiated modified Weibull extension 
distribution 


In this section we introduce the BEMWE distribution using a method similar to that which 
was used by Marshall and Olkin (1967) to define the Marshall Olkin bivariate exponential 
(MOBE) distribution. We start with the joint cumulative function of the proposed bivariate 
distribution and so used it to derive the corresponding joint probability density function. 
Finally The marginal probability density functions and conditional probability density func¬ 
tions of this distribution are also derived. Let X be a random variable has univariate EMWE 
distribution with parameters 7 , a,/3,\> 0 , then the corresponding cumulative distribution 
function (CDF) is given by 


F(x) 


1 — e 


-Aa(e (a,/a) -1) 


X > 0, 


(1) 


and the probability density function (PDF) takes the following form 


/( x) = 7 A/3e (x/a) ' 3 (-) /3 - 1 e- Aa ( e<x/ “ ) ^ 1 ) 
a 


(*/<*)£_ 


1 - e“ Aa(( 


,(x/a)P 


- 1 ) 


-1 7-I 


x > 0 . 


( 2 ) 


2.1 Joint cumulative distribution function 

Suppose that Ui (i = 1,2,3) are three independent random variables such that Ui ~EMWE 
( 7 i ,a, /?, A). Define X\ = max{Ui, U 3 } and X 2 = max{U 2 ,Us}. Then we say that the 
bivariate vector (Xi,Xz) has a bivariate exponentiated modified Weibull extension distri¬ 
bution, with parameters ( 7 i, 7 2 , 73 ) a, P, A) and we denote it by BEMWE( 7l , 72 , 7 3 , a, /?, A). 
The following interpretation can be provided for the BEMWE model. 

Shock model: Assurn thate there exists a three independent sources of shocks. Suppose 
these shocks are affecting a system with two components. It is assumed that the shock 
from source 1 reaches the system and destroys component 1 immediately, the shock from 
source 2 reaches the system and destroys component 2 immediately, while if the shock from 
source 3 hits the system it destroys both the components immediately. Let Ui denote the 
inter-arrival times, between the shocks in source i, i = 1,2,3, which follow the distribution 
EMWE. If X\, X2 denote the survival times of the components, then the bivariate vector 
(Xi, X 2 ) follows the BEMWE model. 

We now study the joint cumulative distribution function of the bivariate random vector 
(X 1 ,X 2 ) in the following lemma. 
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Lemma 2.1. The joint CDF of (.Xl,X 2 ) is 


Fbemwe(xi, X 2 ) 


i - e“ A “ ( 




7i 


1 - e _Aa(e 


72 


1 - e _Aa(e 


73 


where z = min (xi, x 2 ) ■ 


(3) 


proof: Since the joint CDF of the random variables X± and A 2 is defined as 

F(x i,x 2 ) = P (A} < xi, X 2 < x 2 ) 

= P (max{Ui, C/ 3 } < xi, max{[/ 2 , U 3 } < x 2 ) 

= P (Pi < xi,U 2 < x 2 , U 3 < min (xi, x 2 )). 

As the random variables U t (i = 1,2,3) are mutually independent, we directly obtain 

Fbemwe{xi,x 2 ) = P{Ui < xi)P(P 2 < x 2 )P(P 3 < min(xi,x 2 )) 

= F E mwe{x 1 ; 7 i, a, /?, X)F E mwe{x2',12i X )F E mwe(z ; 73 , a, /3, A) (4) 

Substituting from (1) into (4), we obtain (3), which completes the proof of the lemma.2.1. 


2.2 Joint probability density function 

The following theorem gives the joint PDF of the Xi and Af 2 which is the joint PDF of 
BEMWE (7 1 ,7 2 ,7 3 ,a,/3, A). 


Theorem 2.1. If the joint CDF of X\ and X 2 is as in (3), then the joint PDF of X\ and 
X 2 takes the form 


where 

/i(xi,x 2 ) 


/i(xi,x 2 ) if xi < x 2 
IbEMWe{x l,x 2 ) = ^ f2(Xl,X2) if X 2 < Xi 

/o(x,x) if Xi = x 2 = x 


(5) 


/emwe{x 2 ; 7 2 , a , /3, X)fEMWE(xi\‘j 1 + 7 3 , a , /?, A) 

7 2 ( 7 i + 7 s) X 2 f 3 2 e {x2/a)P ( — )/ 3- 1 e - A “( e(X2/ “ > ' 3 - 1 ) 

a 

1 g-Aa(e (x i / “ ),3 -l) i_g-Aa(e <xl/ “ )pf -1) 


1 - e _Aa(e 

71+73-! 


(x 2 /a)/3_ 


7 2 “! 


( 6 ) 


/ 2 (x l,x 2 ) = fEMWE{x\\' 1 i,a,P,X)fEMWE{X 2 \l 2 + l 3 , a iPiX) 
= 7 l (7 2 +7 3 )A 2 /3V a:i/a), 3 (-) / 3 - 1 e- AQ(e(Xl/ “ ),3 - 1) 


a 


xe (®2/a)' 3 j/3-l e -Aa(e <X2/ “^ -1) 

v a 


-Aafe^ 1 ^^ — 1 ) 


1 _ e -Aa(e (x 2 /Q P-l) 


1 — e 

-| 72+73“! 


7i —1 


(7) 
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and 


h{ x -> • 


7s 


7i + 7 2 + 73 


/emwe{x 2 ; 71 +72 + 73: <+ A X) 


7 3 A /3e( x / a )^ (—) ,3_1 e _Aa: ( e(a:/a)/3 - 1 ) 1 — e -Aa(e <x/ “ )/3 -i) 


Ti+72+73- 1 


■ ( 8 ) 


Proof Let us first assume that Xi < X 2 - Then, the expression for fi(xi,X 2 ) can be simply 
obtained by differentiating the joint CDF Fbemwe{x 1 , 2 : 2 ) given in (3) with respect to x\ 
and X 2 • Simillary, we find the expression of f 2 (xi,X 2 ) when X 2 < x\. But f 3 {x, x) can not be 
derived in a similar method. For this reason, we use the following identity to derive f 3 (x, x). 


OO X 2 OO X\ OO 

J J fi(xi,x 2 )dxidx 2 + J J f 2 (xhx 2 )dx 2 dxi+ J f 3 (x,x)d: 


'x = 1 


(9) 


0 0 


0 0 


0 


Let 


OO £2 OOXi 

Ii = J J h( x i,X 2 )dxidx 2 and I 2 = J J f 2 {xi, x 2 )dx 2 dxi 

00 00 

One can find that 

OO 


h= / 72 A/3e (a:2/a),3 ( —) /3_1 e _Aa(e<X2/a>,3_1) 1 - e -\a(e^ /a)fl - 1 ) 

J 2 a 


7i+72+7 3 -! 


dx 2 (10) 


and 


I 2 = [ 7 \p e {xi/a)f> (— 

J a 

0 

Substituting from (10) and (11) into (9) we obtain 

OO 

J f 3 (x, x)dx = 1 - /i - I 2 


1 - e" A “ ( 


,(*l/a) p 


1 7i+7 2 +7 3 — 1 


dx 1 (11) 


Thus, 


OO 

/ (7i+7 2 +7 3 ) A -V 

J « 

0 

oo 

- /7 2 A/3e (a:/a) ' 3 (-) /3 ^ 1 e- Aa(e<x/ “ ) ' 3 ' 1) 

J a 

0 
oo 

J a 


1_ g —Aa(e tx/a ' —1) 


7i+7 2 +7 3 —! 


dx 


l_ e -Aa(e<* / ‘+ 3 -!) 


7i+7 2 +7 3 —1 


dx 


1 _ e -Aa(e<*/“> ,3 -l) 


1 7i+7 2 +7 3 —1 


f 3 (x, x) = i 3 \pe {x/a)f> (±)P- 1 e -W* lm,ar -i) 
a 




I _ g-Aa(e (x/ “ )/3 -l) 


dx. 


1 7i+7 2 +7 3 —! 
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2.3 Marginal probability density functions 

The following theorem gives the marginal probability density functions of X\ and X 2 . 


Theorem 2.2. The marginal probability density functions of Xi ,(i = 1, 2) is given by 


fxAxi) = fEMWE(x i ;'y i +'y 3 ,a,p,X), x t > 0, i = 1,2 


= (7 i + 7s) \pe ^ Xi / a ^-b 


i - e" Aa(e 




-i) 


(7i+7 3 ) — 1 

( 12 ) 


Proof: The marginal cumulative distribution function for X, is 

F(xi ) = P (Xi < Xi) = P ( max{Ui , U 3 } < Xi) = P (Ui < x t , U 3 < Xi ). 

As the random variables Ui (i = 1,2) and U 3 are mutually independent, we directly obtain 


F{xi) 


P (Ui < Xi)P(U 3 < Xi) 

[x _ g-Aa(e (xi/ “ ), 3 -l)j 


1 - e" Aa( 





l _ e -Aa(e (Xi/ “ )/3 -1) 


7i+7 3 


FEMWE^Xi^i + 7 3 , a, /?, A). 


(13) 


From which we readily derive the pdf of Xi , /(#,) = -Jj^F(xi), as in(12). 

2.4 Conditional probability density functions 

The following theorem gives the marginal probability density functions of (Xi,X 2 ). 


Theorem 2.3. The conditional probability density function of Xi given X 3 = Xj , ( i,j = 
1 ,2,* ^ j) is given by 


where 


J x i \x j y a 


fSixA* 


! f Xi \Xji X i I X j) if 0 < Xi < Xj 
fx ) i \x j ( x i\ x i) if 0<Xj<Xi 
fxi\Xj( X i I x i) if X -i = X 3 > 0 


,•) = 


+ 7 3 )^/^ e Xi ' a ( 


(xi/a) f ' | '£ i ' 1 ^-l e -Aa(e (Xi/ “ ), 3 -l) 


1 - e _Aa ( ( 




- 1 ) 


1 7i+7 3 -l 


hi + 73 ) 


1 _ e -Aa(e (x A“)' 3 -l) 


7«+7 3 -l 


,•) = li \p e {xi/a) \—Y- l e- Xa h* i/a)P -V 
a 


1 _ g-AaCe^/^-l) 


1 7,-1 


and 


fxi\x } hi 


73 

7; +73 


1 - e _Aa(e 
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Proof: The proof follows immediately by substituting the joint probability density func¬ 
tion of ( y Xi } X 2 ) given in (6), (7) and (8) and the marginal probability density function of 
given in (12), using the relation 


fXi \XjiXi | Xj) 


fx^XjixuXj) 

fxM) ’ 


1 , 2 ). 


3 Reliability studies 

In this section, we present the joint survival function of (Xi,X 2 ), the CDF of the random 
variable Y = max{Xi, X 2 } and the CDF of the random variable W = min{.Xi, X 2 }. 


3.1 Joint survival function 

In this subsection, we derive the joint survival function of (Xi,X 2 ) in a compact form. 


Theorem 3.1. The joint survival function of (X\. X 2 ) is given by 

Si(xi,x 2 ) if x\ < x 2 
Sx 1 ,Xa(xi,x 2 ) = { s 2 (xi,x 2 ) if x 2 < X\ 

Sq(x,x) if x\ = x 2 = x 


(14) 


where 


Si(xi,x 2 ) = 1- 

x(l- 


l _ e -Aa(e (x2/ “ )3 -l) 


1 7 2 +7 3 


l _ e ~Xa(e ( ' xl/a) - 1 ) 


"I 7i+7 3 


^ _ g-Aa(e^ X2 ' /a ^-l) 


172 


S 2 (X 1 ,X 2 ) = 1- 


1 — e 


—Xoi(e^ Xl ^ ot ^ —1) 


71+73 


X _ e -Aa(e (x 2 /Q °' 3 -l) 


72+73 


X 1- 


X _ g-Aa(e ( * l/ “ )fl - 1 ) 


71 


and 

S 0 {x,x) = 1- 


1 _ e A«(l-e' 


(x/cc)P, 


73 


X _ e -Aa(e (x/ “ ),3 -l) 


1 7i 


1 _ e -Aa(e^“)"-l) 


7 72 


X _ e -Aa(e^“+-l) 


1 71+72 


Proof: The joint survival function of (Xi, X 2 ) can be obtained from the following relation 
Sx 1 ,x 2 (xl,x 2 ) = 1 - F Xl (xi) - F X a(x 2 ) + F Xlt x a (xi,X 2 ). (15) 
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Substituting from (3) and (13) in (15), we get 


S Xl ,x 2 (xi,x 2 ) = 1 - 


1 - e“ Aa(e 


(xi/a )/ 3 


-1) 


1 7i+7 3 


1 - e“ Aa(e 


(x 2 /a)/3 


-1) 


-| 72+73 



7i 

J _ g -Aa(e (a:2/ “ ) ' 3 -l) 

72 

l _ g-Aa(e ( * / “ )/3 -1) 

. 


. 


. 


16 ) 


where 2 = min (xi, X 2 ) ■ From (16) we can be obtained simply the expressions of S\(x 1 , £ 2 ), 
S 2 (xi,x 2 ) and So(xi,x 2 ) for x\ < X 2 , X 2 < X\ and X\ = X 2 = x respectively, which 
completes the proof. 


Comment 3.1. Basu (1971) defined the bivariate failure rate function h(x 1 , 2 : 2 ) for the 
random vector (Xi,X 2 ) as the following relation 


hx 1 ,X 2 {xi,X 2 ) 


fx 1 ,x 2 (xi,x 2 ) 
Sx 1: x 2 (Xl,X 2 ) 


(17) 


We can obtained the bivariate failure rate function h(x 1 , X 2 ) for the random vector ( X \, X 2 ) 
by substituting from (5) and (14) in (17). 


Lemma 3.1. The CDF of the random variable Y = max{X\, X 2 } is given as 


Fy(y) 


Y Xa(e^ v ^ a ^ — 1) 


7i+7 2 +7 3 


(18) 


Proof: Since 

F y (v) = P(Y<y) = P(max{X 1 ,X 2 }<y) = P(X 1 <y,X 2 <y) 

= P ( max{Ui , U 3 } < y, max{U 2 , U 3 } <y) = P{Ui< y, U 2 <y,U 3 <y), 

where the random variables Ui (i = 1,2,3) are mutually independent, we directly obtain 

Fy(y) = P(U 1 <y)P(U 2 <y)P(U 3 <y) 

= F E MWE{y,'yi,ot,p,X)FEMWE(y]'y2^ a ^^)F E MWE(y,'y 3 ,ot,p,X). (19) 

Substituting from (1) in (19), we get (18) which completes the proof of the lemma.3.1. 


Comment 3.2. From lemma 3.1. we can say that, if X\ and X 2 are independent EMWE 
random variables then max{Xi,X 2 } is also EMWE random variable. 


Lemma 3.2. The CDF of the random variable W = rriin{Wi, W} is given as 


F w {w) 


l _ g-Aa(e ( “’ / “ )/3 



1 - e“ Aa( 






(w/ a p 


- 1 ) 


7i+7 2 +7 3 


( 20 ) 
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Proof: Since 

Fw{w) = P (W < w) = P (min{X 1 ,X 2 } < w) = 1 — P (min{X 1) X 2 } > w) 

= 1 — P > w, X 2 > w) = 1 — S(w, w) (21) 

Substituting from (14) in (21), we get 

Fwiw) = F Xl (w) + F X2 (w) - F x u x 2 {w, w). (22) 

Substituting from (3) and (13) in (22), we get (20) which completes the proof of the lemma 
3.2. 


4 The marginal expectation 

In this section, we derive the marginal expectation of Xi (i = 1, 2). The following theorem 
gives the rth moments of Xi (i = 1, 2) as infinite series expansion. 


Theorem 3.1. The rth moment of Xi (i = 1, 2) is given by: 


OO OO 




av 


j —0 k —0 


( 7 i + 73) - A (-l) J+k X k a k+l3+r U + l) k 
j ) (k +1 y+ r k\ 


e \ a (j+ l) r (!_ + 

(23) 


Proof: We will start with the known definition of the rth moment of the random variables 
Xi with pdf f(xi) given by 

OO 

E(XT ) = J x r i fx i {x i )dxi. 

0 

Substituting for from (12), we get 


OO 

E ( X r ) = (7 i + 7 3 )-A J x r+ 0 -l e ( Xi/a )e e -\ a (e^/^- 1 ) 


l _ e ~\a(e < - Xi/a ' ,f> -1) 


1 ( 7 .+ 73)- 1 


dXi- 


(24) 

Since 0 < e~ Xa ^ e< ' ,X ' < 1 for x > 0, then by using the binomial series expansion of 


1 - e~ Xa (' 


( 7 «+ 73>- 1 . 


- 1 ) 


given by 


l _ - 1 ) 


= “ A7l + 73) _ 1 
3=0 


J 


(-l) J e 


,3 _-jAa(e <x 7“) _i) 


(25) 


Substituting from (25) into (24), we get 

E(XV ) = Ai + 7s) A f A* + ^3) - A ( -_ 1 y e \a(j+l) f x r+h-l e {x i / a ^ e -\ a U+l)e < - x P a '> P 

a j =0 \ 3 / J 


0 



Using the series expansion of e 


-Aa(j + l)e <x i /Q °' 3 


, one gets 


mi) = 


hi + 7 3 ) 

aP~ l 


EE 


j =0 k =0 
00 

x fxh^e^+^/^dxi. 


hi + 73) - A (-1) J+A A fc q fc (j + l) fc Aa(.7 +1 ) 
j ) k\ 


Let y = (k + \)(xi/&Y i n the above integral, then we can get 

hi + 7s) - A hh 1+k A k a k +P+ r (j + l) fe 


OO OO 


T?(vr \ _ hi + 73) A \ ' \ ' 

E ( Xi ' a P-l 


j —0 k—0 


{k + 1 )P+ r k\ 


OO 

e Aa(i+l) J e y dy 


(26) 


OO 

Since, T(z) = x z J e xt t z ~ l dt , z > 0, x > 0, then 


OO 

J yTe v dy = Y(k + 1). 


(27) 


Substituting from (27) into (26), we get (23). This completes the proof. 


5 Maximum liklihood estimators 

In this section, we use the method of maximum likelihood to estimate the unknown param¬ 
eters of the BEMWE distribution. Consider constant values to the parameters a, [3 and A 
so, we want to estimate the other parameters an d 73- Suppose that we have a sample 

of size n ,of the form {(a:ii,X 2 i), (£ 12 , X 22 ),•••, (xi n ,x 2n )} from BEMWE distribution. We 
use the following notation 

h = {xu<X 2 i}, h = {xu > x 2i }, h = {xu = x 2i = Xij, / = JiU/ 2 U/ 3 , 
|/i|=rii, \I 2 \ = n 2 , |/ 3 | = n 3 , and ni + n 2 + n 3 = n. 

Based on the observations, the likelihood function of the sample of size n given by: 

ni ri2 n 3 

j(7l> 72) 73> «, h A) = fl(xii,X 2 i) P f 2 {xii, X 2 i) P f 3 (Xi, x). 

i— 1 2=1 2=1 

The log-likelihood function can be written as 


L( 7 i> 7 2 > 73> h A) = n 1 In (j 2 {-y 1 + 73) \ 2 (3 2 ) + E(“) /? “ Xa ^{e {xii/a)P - I) 
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'«■ I '‘■I 0 "-l 

-AaE( e(W< M - 1) + (72-1) E In(l - g-W 2 * 7 * 0 ->) + (0 - 1) Eln(^) 

2 = 1 2=1 2=1 

ni ni 

+ (7i+ 73- 1) E ln(l - e -Aa(e(-H /“)' 3 -1)) + + n 2 ln ( 7l ( 72 + 73 ) A 2 /3 2 ) 

2=1 2=1 ^ 

222 n 2 212 212 

+ f (-f - Aa Vte' 1 '*/ 11 )' _ 1) (-f -XaY(e {x2i/a)f> -1) 

2=1 2=1 2=1 2=1 

n 2 0 21 2 

+ ( 7l -l) £ ln(l - e-W* 1 ’ 1 *'* -V) + (72+73-1) E Ml “ e~^ 2i/a) ^) 

2=1 2=1 
212 213 2I3 

+ (/?-!) E M^r) + n 3 ln (7s A /3) - E ( e( " i/a)/3 “ + E (-)' 3 

1 ^ -1 1 ^ 

2=1 2=1 2=1 

223 213 ^ 

+ (/3 - !) E M 77 ) + (7i+7 2 + 7 3 _ 1 ) E Ml _ e _A “ (e<a! ‘ /a) _1) ). (28) 


i—1 


i=l 


Computing the first partial derivatives of (28) with respect to 71,72 and 73 and setting the 
results equal zeros, we get the likelihood equations as in the following form 


dL 

dli 


n 1 


7i +73 


21 1 a 21 2 

+ E ln(l - e- A “+ <Xli/ “ ) _1 )) + — + Eln(l - e - A “< e< ' li/o> ~V) 

2=1 «-1 


2=1 


+ E ln (! - 1 


(29) 


dL 

dl2 


2=1 


n l , i._/i 0 -Aa(e (x 2i/“>' 3 ~l), 




ri 2 


72+73 


EM 1 - 


— Aa( 




-!)1 


«3 


EMI- 


Aa(i 


.(x./a)/ 3 , 


(30) 


i=l 


and 


dL n 1 

^73 7i + 73 


EMI- 


, —Aa( 


_i 




72+73 


EMi-' 


Aa(' 


>( x 2i/ a ) P _ 


+ ^ + 

7s 


E mi - 


, —Aa(e' 


( x i/<*) p 


- 1 )! 




(31) 


i=1 


To get the MLEs of the parameters 7l , i 2 and 73 , we have to solve the above system of 
three non-linear equations. The solution of equations (29), (30) and (31) are not easy to 
solve, so numerical technique is needed to get the MLEs. 


5.1 Asymptotic confidence bounds 

In this subsection we consider the approximate confidence intervals of the parameters 7l , 7 2 
and 73 by using variance covariance matrix 7 q see Lawless (2003), where I ^ 1 is the inverse 
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of the observed information matrix 



( d 2 L 

3rf 

d 2 L 
d7i372 

a 2 L \ 
37i37 3 

-1 

( Var{ 77 ) 

Cov{ 7 i,7 2 ) 

Cow( 7 i, 7 3 )\ 

lo 1 = - 

d 2 L 
dl 2 37 1 

d 2 L 

d 2 L 

372373 

= 

Cov{ 7 2 , 7 i) 

Va,r( 7 2 ) 

C°v{ 72 , 73 ) 


d 2 L 

\d7 3 d7i 

a 2 l 

373 373 

d 2 L , 

37s / 


\Ccw(73,7i) 

Cov (73,72) 

V(ir{ 73) ) 


The derivatives in I 0 1 are given as follows 


d 2 L 

ni n 2 

92 L -0 

d 2 L m 

Wi 

(7i+73) 2 7i 2 ’ 

971^72 ’ 

^71^73 (7i+73) 2 


d 2 L 

ni 

n 2 

d 2 L n 2 

W 2 

72 

(72 +73) 2 ’ 

^72^73 (72 + 7 3 ) 2 


and 

d 2 L Hi 712 7X3 

dl l (7i + 73) 2 (7 2 + 7s ) 2 7 2 ' 

We can derive the (1 — 5)100% confidence intervals of the parameters 7 l5 7 2 and 73 by using 
variance covariance matrix as in the following forms 



where Z s is the upper (|)th percentile of the standard normal distribution. 


6 Data analysis 

In this section we present the analysis of a data set and we consider a constant value to the 
parameters a, 13 and A which take the values 0.1, 0.3 and 0.05 respectively. The data set 
has been represent the American Football (National Football League) League data and they 
are obtained from the matches played on three consecutive weekends in 1986. The data 
were first published in ‘Washington Post’ and they are also available in Csorgo and Welsh 
(1989). It is a bivariate data set, and the variables Xi and X 2 are as follows: X\ represents 
the ‘game time’ to the first points scored by kicking the ball between goal posts, and X 2 
represents the ‘game time’ to the first points scored by moving the ball into the end zone. 
These times are of interest to a casual spectator who wants to know how long one has to 
wait to watch a touchdown or to a spectator who is interested only at the beginning stages 
of a game. 

The data (scoring times in minutes and seconds) are represented in Table 1. Here also 
all the data points are divided by 100 just for computational purposes. The variables have 
the following structure: (i) X\ < X 2 means that the first score is a field goal, (ii) X\ > X 2 , 
means the first score is an unconverted touchdown or safety, (iii) X\ = X 2 means the first 
score is a converted touchdown. 
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Xi 

*2 

Xi 

a 2 

X x 

a 2 

Xx 


2.05 

3.98 

8.53 

14.57 

2.90 

2.90 

1.38 

1.38 

9.05 

9.05 

31.13 

49.88 

7.02 

7.02 

10.53 

10.53 

0.85 

0.85 

14.58 

20.57 

6.42 

6.42 

12.13 

12.13 

3.43 

3.43 

5.78 

25.98 

8.98 

8.98 

14.58 

14.58 

7.78 

7.78 

13.80 

49.75 

10.15 

10.15 

11.82 

11.82 

10.57 

14.28 

7.25 

7.25 

8.87 

8.87 

5.52 

11.27 

7.05 

7.05 

4.25 

4.25 

10.40 

10.25 

19.65 

10.70 

2.58 

2.58 

1.65 

1.65 

2.98 

2.98 

17.83 

17.83 

7.23 

9.68 

6.42 

15.08 

3.88 

6.43 

10.85 

38.07 

6.85 

34.58 

4.22 

9.48 

0.75 

0.75 



32.45 

42.35 

15.53 

15.53 

11.63 

17.37 




Table (1). American Football League (NFL) data 


From this data, we find that the values of the unknown parameters 7 ), q 2 and 7^3 are 
0.0416, 0.253 and 0.52 respectively and the log-likelihood equals (-250.28 ). By substituting 
the MLE of unknown parameters in (32), we get estimation of the variance covariance matrix 
as 

/ 0.000842 0.00000395 -0.000299 

/o ' 1 = 0.00000395 0.00394 -0.0000939 

\-0.000299 -0.000094 0.00711 

The 95% confidence intervals of 7 ), 7 2 and 73 are (0,0.098), (0.130,0.376) and (0.355, 0.685) 
respectively. 



7 Conclusions 

In this paper we have introduced the bivariate exponentiated modified Weibull extension 
distribution whose marginals are exponentiated modified Weibull extension distributions. 
We discussed some statistical properties of the new bivariate distribution. Maximum likeli¬ 
hood estimates of the new distribution are discussed and we provided the observed Fisher 
information matrix. One real data sets are analyzed using the new distribution. 
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